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ABSTRACT 

We report the results of high spatial and spectral resolution integral-field spectroscopy of the 
central ~ 3 x 3 arcsec 2 of the active galaxy NGC 1275 (Perseus A), based on observations 
with the Near-infrared Integral Field Spectrograph (NIFS) and the ALTAIR adaptive-optics 
system on the Gemini North telescope. The circum-nuclear disc in the inner R ~ 50 pc of 
NGC 1275 is seen in both the H2 and [Fe II] lines. The disc is interpreted as the outer part 
of a collisionally-excited turbulent accretion disc. The kinematic major axis of the disc at a 
position angle of 68° is oriented perpendicular to the radio jet. A streamer-like feature to the 
south-west of the disc, detected in H2 but not in [Fe 11], is discussed as one of possibly several 
molecular streamers, presumably falling into the nuclear region. Indications of an ionization 
structure within the disc are deduced from the He I and Br7 emission lines, which may par- 
tially originate from the inner portions of the accretion disc. The kinematics of these two lines 
agrees with the signature of the circum-nuclear disc, but both lines display a larger central ve- 
locity dispersion than the H2 line. The rovibrational H2 transitions from the core of NGC 1275 
are indicative of thermal excitation caused by shocks and agree with excitation temperatures 
of ~ 1360 and ~ 4290 K for the lower- and higher-energy H2 transitions, respectively. The 
data suggest X-ray heating as the dominant excitation mechanism of [Fe II] emission in the 
core, while fast shocks are a possible alternative. The [Fe II] lines indicate an electron density 
of ~ 4000 cm~ 3 . The H2 disc is modelled using simulated NIFS data cubes of H2 emission 
from inclined discs in Keplerian rotation around a central mass. Assuming a disc inclination 
of 45° ± 10°, the best-fitting models imply a central mass of (8 ±2) x 10 8 M Q . Taken as a 
black-hole mass estimate, this value is larger than previous estimates for the black-hole mass 
in NGC 1275, but is in agreement with the M-a relation within the rms scatter. However, 
the molecular gas mass in the core region is tentatively estimated to be non-negligible, which 
suggests that the central mass may rather represent an upper limit for the black-hole mass. In 
comparison to other H2-luminous radio galaxies, we discuss the relative role of jet feedback 
and accretion in driving shocks and turbulence in the molecular gas component. 

Key words: galaxies: general - galaxies: active - galaxies: nuclei - galaxies: individual 
(NGC 1275) 



1 INTRODUCTION 

NGC 1275 (also commonly known as Perseus A, 3C 84, 4C+41.07, 
or IRAS 03164+4119) is the central galaxy in the Perseus clus- 
ter, Abell 426. The cluster Abell 426 shows two distinct ve- 
locity components: NGC 1275 (z = 0.017559) belongs to the 
low-velocity component, which is associated with velocities of 
5200 kms -1 . The high-velocity component, showing velocities 
of 8200 kms- 1 is located in front of NGC 1275 to the north 
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and north-west of the nucleus dRubin etal.1ll977l) . The two ve- 
locity components ha ve been interpreted as indicating an ongoing 
galaxy cluster merger Jjjnger et al. Il990l). However, X-ray absorp- 
tion studies bv lGillmon et al. I J2004l) and lSanders & Fabian! ( 2007) 
indica te that the separation between both velocity components is 
large. iGillmon et al. 1 J2004I) suggest that the high-velocity system 
is not yet merging with NGC 1275 but is still in the process of 
falling into the cluster. 

The number of classifications for NGC 1275 in the literature - 
Fanaroff-Riley 1 radio galaxy (F-R I), BL Lac object, post-merger 
galaxy, cooling-flow galaxy, peculiar Seyfert 1, and LINER - re- 
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fleet the fact that NGC 1275 exhibits an extraordinarily rich set 
of interstellar physical phenomena. It is a star-forming early-type 
cD galaxy, h arbouring an active nucleus which is a poi nt source of 
hard X-rays dKawara & Taniguchilll993l : |Prietolll996t) and a very 
strong radio source (hence the name Perseus A). The ra dio emission 
on large-scales and parsec-scales shows a pair of jets jPedlar et al. I 
Il990l : Walk er. Romnev & Benson II 19940. T he jets are oriented at a 



position angle of 160° dPedlar et al. II 1990]). Using VLBI measure- 
ments of the parsec-scale jet. Iwalker. Romnev & Benson 1 Jl994h 
inferred a likely range of angles to the line of sight of 6 = 
30° — 55°. The southern jet is the one approaching us. 

NGC 1275 is associated with strong thermal X-ray emis- 
sion, which extends over several arc minutes around the active nu- 
cleus and sho ws a complex bubble-like morphology I Fabia net al. I 
l2000l [2003a), as well as with an extended structure of Ha fila- 
ments ( Min kowski II 1 9571 : iLvndsll 1 970p . These filaments of ionized 
gas ha ve been studied in detail by IConselice. Gallagher & Wvse I 
2001) . They are likely to be stabilized by magnetic fields 



(Fabia n et al. 11200^) and display L I NER-like spectra w hich have 
been analysed bv lSabra et al. ll200oh . lSabra et al. I feood) examined 
a variety of excitation mechanisms - including shocks, ionization 
by an AGN, and stellar photoionization - without finding a fully 
adequate explanation for the observed spectra. Collisional heating 
by ionizing particles has recently been investigated in greater de- 
tail as a po ssible explanation for the emission char acteristics of the 
filaments l lFerland et al. 1120091 : iFabian et al. 11201 ll) . 

NGC 1275 was originally believed to represent the prototypi- 
cal cooling-flow galaxy, in which t he Ha filament s have formed by 
cooling of the intracluster medium jFabian II 19941) . But the correla- 
tion between the X-ray and Ha morphology suggests that the fila- 
ments rather consist of cold gas which is uplifted b e hind the radio 
bubble s dri ven by the rad i o jet jFabian et alj|2003l) . TSalome et al.l 
(2006) and Salo me et al. I J201 ll) have shown that there is also a 
tight correlation between the extended CO line emission, the Ha 
filaments, t he extended X-ray emission, an d the large radio lobes 
mapped by Pedl ar et al. I dl990h . Likewise. iLim et all 12012) found 
a close correspondence between the Ha filaments and the extended 
H2 1-0 S(l) emission. A fraction of the molecular gas and much of 
the thermal X-ray emission are confined in shells around bubbles 
of non-thermal synchrotron-emitting plasma. This is clear evidence 
of the importance of the radio lobes in compressing and cooling the 
ambient galactic medium. This ambient galactic medium can be de- 
rived from either shock-cooled the rmal X-ray plasma or from the 
debri s of a recent galactic merger dConselice. Gallagher & Wvse I 
120011) . 

The total mass of molecular gas in NGC 1275 is large. 
Based on CO(2-l) o bservations of the nuclear and extended gas, 
ISalome et alj d2006t) measured a total mass of cold H2 of 4 x 
10 10 M Q , using the standard CO to molecular hydrogen conver- 
sion factor. About 3 x 10 9 Mq of this molecular gas are contained 
in the inner R ~ 3 kpc, a mass that is comparable with the to- 
tal mass of H2 in the Milky Way. Most of the emission arises 
in an extended eas t-west structure passing through the nucleus. 
ISalomeetal]d2006h found that the CO emission in the central 8 kpc 
diameter region does not show any evidence of large-scale rotation. 
However, the CO line profile in the central 4 kpc diameter region is 
doubl e-peaked, which ma y indicate inflows, outflows, or a rotating 
disc dSalomeetaill201ll) . 

The source of the molecular hydrogen in the nuclear region of 
NGC 1275 is thought to be turbulent infalling gas. Whether this gas 
is carried into the ce ntral regions by a g alactic merger or a cooling 
flow is unclear (e.g. ISalome et al. 2006). Interferometric observa- 



tions of the CO emission in the central 10 kpc bv lLim etal.1d2008l) 
revealed that most of the gas is associated with radial filaments, 
which are oriented in the east-west direction (i.e. perpendicular to 
the jet axis), located betwe en X-ray cavities, and consistent with 
radially infalling gas flows. ILim et al~l d2008l) argued that the ab- 
sence of any large-scale rotation makes a merger-origin of the gas 
less likely. 

The near-infrared spec trum of the nuclear region was inves- 
tigated by iRudv et ail d 19931) who found that the [Fe II] lines are 
unusually strong and narrow, indicating that they arise from a 
region with kinematics different from the nucleus. It was sug- 
gested that these features are caused by shocks. Shocks were 
also suggested as the ca use of the H2 1-0 S(l) line emis- 
sion, first investigated by iFischer et al.l dl987l) . These molecu- 
lar sh ocks were discussed as possible triggers of starburst ac- 
tivity. iKrabbe et al. I d2000h. iRodriguez-Ardila. Riffel & Pastorizal 
(2005), and Wilman, Edge & Johnstone I d2005l) determined the ex- 
citation temperatures of the near-infrared H2 lines in the circum- 
nuclear region. They found the emission of this ~ 1300 — 4000 K 
H2 component to be consistent with thermal excitation from X-ray 
or shock heating. Because th e brightest H2 emiss ion is extended 
over only the central ~300 pc jKrabbe et al. I d2000l) concluded that 
it is more likely to be directly associated wi th the AGN rather than 
being located in a large-scale cooling flow. iDonahue et al. I j200Ct) 
agreed with this conclusion and further suggested that the H2 emis- 
sion is too bright to be part of a cooling flow. In addition, their HST 
images show a faint extension of H2 1-0 S(l) emission from the 
nucleus to the southwest as well as a faint feature to the east of the 
nucleus. 

Wilman, Edge & Johnstone 1 d2005h inferred the properties of 
the near-infrared H2 in the centre of NGC 1275, using observa- 
tions of the central 6.5 x 3.4 arcsec 2 with the imaging spectro- 
graph UIST (UKIRT Imaging Spectrometer) on UKIRT (United 
Kingdom Infrared Telescope) taken under seeing-limited condi- 
tions of 0.3-0.4 arcsec. They suggested that the H2 emission at 
around 50 pc from the nucleus arises in a dense clumpy disc, rotat- 
ing in a plane perpendicular to the radio jet. They concluded that 
thermal excitation by the nuclear radiation field is the dominant 
excitation mechanism of the H2 emission. Based on the discon- 
tinuity in H2 velocity across the nucleus of around 240 km s~ , 
they estimated a mass of 3.4 x 10 8 M© for the black hole at 
the nucleus, assuming a disc inclination of i = 45°. The nu- 
clear [Fe II] 1.644 [im emission was found to be spatially resolved, 
w hile the nuclear Paq emission rema ined unresolved. According 
to Iwilman. Edge & Johnstone 1 d2005l) . the nuclear disc of hot H2 
would likely be part of the larger-scale H2 and CO structure seen 
in NGC 1275. This, together with the more recent observations of 
molecular gas in the extended filam ents on several 10 kpc scales 
( Salo me et al ,l2006l : lLim et alj2012h . suggests that the fuelling pro- 
cess of the AGN in NGC 1275 by interstellar material is being ob- 
served over a large range of scales. 

In this paper, we analyse the inner circum-nuclear disc of 
NGC 1275 using o bservations with NIFS (Near-infrared Integral- 
Field Spectrograph; iMcGregor et al. l2003h in conjunction with the 
adaptive-optics system, ALTAIR, on the Gemini North telescope. 
The NIFS data permit a high-spatial-resolution study of the kine- 
matics of the molecular and ionized gas via the H2 1-0 S(l), [Fell] 
1.644 \im, He I 2.058 |^m, and Br7 2.166 [im lines as well as an ex- 
amination of the excitation mechanisms of H2 and [Fe II], Through- 
out the paper we assume Ho = 71 km s _1 Mpc -1 , Qm = 0.3, and 
!!,„ = 0.7 , which results in a spatial scale of 358 pc arcsec -1 for 
NGC 1275 dWrightll2006l) . 
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2 OBSERVATIONS 

The Gemini Near-infrared Integral-Field Spectrograph (NIFS) and 
ALTAIR adaptive-optics system on the 8.1m Gemini North tele- 
scope on Mauna Kea in Hawaii were used to obtain ii-band and 
it -band spectra of NGC 1 275 . NIFS has 29 slitlets, each 0. 1 arcsec 
wide, and has a resolution of 0.04 arcsec per spaxel along each slit- 
let. The field of view is 3.0 x 3.0 arcsec 2 , and it delivers a spectral 
resolving power of ~ 5300, matched to two pixels on the detector, 
at both H and it. This corresponds to a velocity resolution of ~ 
60kms _1 . 

The A"-band observations of NGC 1275 were obtained on 
2005 November 1 1 during the commissioning of NIFS. The obser- 
vations were performed using the ALTAIR natural guide star sys- 
tem with the nucleus of NGC 1 275 as the reference for the adaptive- 
optics correction. The nearby star U1275.02199239 was used as 
the on-instrument wavefront sensor (OIWFS) guide star for slow 
flexure correction. The instrument was oriented at a position angle 
of 90°, so that the slitlets were aligned E-W. The observation con- 
sisted of eight exposures of 600 s duration centred on NGC 1275, 
interleaved with eight 600 s exposures of an adjacent sky region. 
None of these exposures was spatially dithered. The A0V HIPPAR- 
COS star HIP 18769 was observed after the NGC 1275 observa- 
tion to provide correction of telluric absorption and flux calibration 
from its 2MASS magnitudes. 

The ii-band observations of NGC 1275 were obtained during 
NIFS Science Verification observations (program GN-2006A-SV- 
124) on 2006 February 1 and 2. Again, ALTAIR was used in natu- 
ral guide star mode and the instrument was set to a position angle 
of 90° . Two object-sky exposure pairs were obtained with 900 s ex- 
posure times in modest seeing, one pair on each night. HIP 18769 
was also observed after NGC 1275 on each night. 



3 DATA REDUCTION AND VISUALIZATION 

The d ata were reduced in a standard way dStorchi-Bergmann et al. I 
12009b using a variant of the Gemini IRAF package. The object and 
sky exposures were combined first. Then a dark exposure was sub- 
tracted from both the combined-object and combined-sky frames. 
The combined-sky frame was scaled to match the combined-object 
frame and subtracted. The sky-subtracted object frame was then 
flatfielded, bad pixels were interpolated, and the 29 individual two- 
dimensional slitlet spectra were extracted from this image. Ar- 
gon arc spectra were used to define the wavelength calibration 
and Hatfield exposures through a Ronchi grating were used to de- 
fine the spatial calibration along slitlets. Each two-dimensional 
slitlet spectrum was then geometrically transformed onto a rec- 
tilinear spatio-spectral coordinate grid and the resulting 29 two- 
dimensional spectra combined into a three-dimensional data cube. 
A one-dimensional spectrum was extracted from a 0.5 arcsec - 
diameter circular aperture defined within the standard star data 
cube, a correction for telluric absorption was defined, and this one- 
dimensional correction spectrum was applied to all spaxels within 
the object and standard star data cubes. Similarly, a flux-calibration 
spectrum was extracted from the telluric-corrected standard star 
data cube using a 1.5 arcsec-diameter circular aperture with the 
level set by the 2MASS H or it band magnitude and the black- 
body shape defined by its 2MASS J — K colour temperature of 
1 1 000 K. The final data cubes obtained at H and it are in the 
form of multi-extension FITS files with 2040 wavelength pixels 
and 29 x 69 spatial pixels (spaxels). 
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Figure 3. PSF fit (solid line) to the it-band profile of the standard star 
HIP 18769 (plus signs) and a comparison with the profiles of the it-band 
continuum (filled squares) and the integrated Br7 emission (open squares) 
of NGC 1275. The Br7 emission is integrated over a velocity interval from 
—500 km s" 1 to +500 km s — 1 . All profiles are centred on the emission 
centroids, respectively, and normalized to their peak values. Positive and 
negative radii show the profile obtained from an azimuthal average over 
two opposite halves of the emission peak. The PSF fit is based on a Moffat 
profile, as explained in the text. 

The resulting H- and it-band spectra of the core of NGC 
1275, integrated over a 1.0 arcsec-diameter circular aperture, are 
shown in Figs Q] and [2] The it-band spectrum of the core of 
NGC 1275 is known to be dominated by H2 emi ssion lines 
dKrabbe et al. l200d ; lwiiman. Edge & Johnstone 1 2005). Our spec- 
tra confirm a large number of H2 lines, mostly in the it-band spec- 
trum and some weaker lines in the ii-band spectrum. The ii-band 
is dominated by the strong [Fe II] 1.644 [im emission. Other weaker 
[Fe II] emission lines are also apparent. Furthermore, weaker lines 
of the Brackett series and of He I are seen in the H- and it-band. 
The flux measurements for the emission lines detected in the NIFS 
spectra are listed in TableQ] 

A knowledge of the it-band point-spread function (PSF) 
is required for analysing the surface-brightness profiles and the 
kinematics of the emission-line gas in the circum-nuclear region 
of NGC 1275. We use the spatial profile of the standard star, 
HIP 18769, as our PSF reference. Obtaining the PSF from sepa- 
rate standard-star observations involves uncertainties, since the ob- 
serving conditions may have changed between these and the ob- 
ject observations. The PSF could be measured more directly using 
an unresolved emission component in the it-band data cube for 
NGC 1275. Candidates for this are the it-band continuum emis- 
sion or hydrogen emission from the AGN broad-line region. How- 
ever, Fig. [3] shows that the profile of the it-band continuum of 
NGC 1275 (filled squares) is extended compared to the profile of 
the standard star. The same is true for the Br7 2. 166 [im emission of 
NGC 1275 (open squares), which is integrated over a velocity inter- 
val from —500 km s~ x to +500 km s" 1 after continuum subtrac- 
tion. The fact that the Br7 2.166 \im emission is spatially extended 
suggests that the integrated line includes emission from outside the 
AGN broad-line region. In order to isolate any emission from the 
AGN broad-line region, it would be necessary to decompose the 
line profile, which is not practical at the given signal-to-noise ratio. 

The radial profile of HIP 18769 and the corresponding PSF fit 
are shown in Fig.[3] The PSF fit is based on a single Moffat function 
of the form 



© RAS, MNRAS 000,[TJfT9] 



4 J. Scharwachter et al. 



Table 1. Emission-line fluxes for integrated spectra centred on the nucleus of NGC 1275. Column 1: Name of transition; Column 2: Vacuum wavelength of the 
line; Columns 3 and 4: Flux and Gaussian FWHM of the lines in the spectrum integrated over a 1 arcsec-diameter aperture centred on the nucleus; Column 5 
and 6: Flux and Gaussian FWHM of the lines in the spectrum integrated over a 0.5 arcsec-diameter aperture centred on the nucleus. Column 7: Comment 
regarding unlisted values. The errors of the flux and FWHM measurements are both estimated to be 20 per cent. 
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where Co is the amplitude, R is the radius, and FWHM is the full 
width at half maximum. Fitting this function to the standard star 
profile results in a FWHM of 0.103 arcsec (36 pc), corresponding 
to one spatial pixel across the slitlets. 



4 EMISSION-LINE MORPHOLOGY 

The circum-nuclear regions of the integrated flux maps in H2 1-0 
S(l) and [Fe II] 1.644 \im show a steeply increasing surface bright- 
ness towards the core of the galaxy (Fig.|4^ and c). In the analysis 
presented here, this emission is interpreted as tracing the outer parts 
of an inclined accretion disc around the central active nucleus. The 
classical inner accretion disc is unresolved at our spatial resolution 
of 36 pc. 

The integrated H 2 emission is clearly perturbed towards the 
south-west. Fig.|4p shows an unsharp mask of the integrated H 2 
image of Fig. [4^, which is the result of subtracting a smoothed 
copy of Fig.|4^ from Fig.|4^. The smoothing is performed by con- 
volving the image with a Gaussian function, using a FWHM of 
3 pixels. Fig. |4p reveals a feature extending from the nucleus to 
the west and south, which is like ly to be a part of the s outhwest- 
ern emission feature reported bv lDonahueetal.l bOOfj) . It is un- 
certain whether this is a stream of gas fa lling into the ce ntre or 
being ejected from the centre. However, iLim et al~l {2008) found 
signs of generally infalling motions in the molecular gas, CO(2-l), 



in the inner R ~ 8 kpc on kiloparsec scales. Since it is likely that 
these streams of gas continue to smaller scales, we interpret the 
southwestern feature in our NIFS data as an infalling stream of 
gas settling into the overall accretion disc at small radii. As we 
will see in Sections 151 and [8^21 the streamer is associated with red- 
shifted velocities and apparently joins the H 2 disc on a retrograde 
orbit. The streamer is not detected in [Fell] 1.644 |j.m emission. 
This is obvious from Fig. |4ji, which shows the unsharp mask cor- 
responding to the integrated [Fell] 1.644 \im emission. The lack 
of [Fe II] detection in the streamer suggests that the streamer is ei- 
ther predominantly molecular in nature, or not exposed to sufficient 
internal high-velocity shocks or sufficient X-ray heating from the 
nuclear continuum source, both of which are possible [Fe II] exci- 
tation mechanisms. 



5 EMISSION-LINE KINEMATICS 

In agreement with the morphological evidence for a circum-nuclear 
disc from Section|4]and the previous i ndications for a rotating disc 
based on the lower-resolution data bv lWilman. Edge & Johnstone I 
d2005h . the H 2 1-0 S(l) 2.122 11m and [Fell] 1.644 ixm emis- 
sion lines show kinematics indicative of disc rotation in the central 
R ~ 0.15 arcsec (~ 50 pc) of NGC 1275. The kinematic maps, 
based on a single-Gaussian profile fit to the H2 1-0 S(l) 2.122 \im 
and [Fell] 1.644 \im lines, are shown in Figs [5] and [6] The four 
panels of each figure show the central velocity (upper left panel), 
the velocity dispersion (upper right panel), the continuum flux den- 
sity (lower left panel), and the intensity of the Gaussian fit to the 
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Figure 1. //-band spectrum of the core of NGC 1275 integrated over a 1.0 arcsec-diameter circular aperture. The wavelength scale is shown as observed for 
NGC 1275 at z = 0.017559. Note that the [Fe II] lines are broader than the H 2 1-0 S(7) transition. 



line profile (lower right panel). The scenario of disc rotation is sup- 
ported by the fact that the redshifted and blueshifted emission of 
H 2 1-0 S(l) 2.122 [im and [Fell] 1.644 is symmetrical with 
respect to the peak in continuum intensity (indicated by the cross 
in Figs. [5] and [6}. In addition, the continuum peak coincides with 
the line intensity peak and the peak in velocity dispersion of both 
lines. 

A jet-driven origin of the H2 and [Fell] kinematics in the 
circum-nuclear region can be ruled out, because the major kine- 
matic axis of the proposed disc is perpendicular to the radio jets. 
The major kinematic axis is found at a positi on angle of PA=68° , 
while the radio jets are located at PA=160° jPedlar et ai7lll990h . 
The fact that the two position angles are perpendicular to each 
other strongly suggests that the observed rotating disc represents 
the outer parts of an accretion disc and that the jets are ejected 
along the polar direction. In this case, a disc inclination of about 
45° can be inferred from the known angle of the radio jets wit h 
respect to the line of sight (cf. lWilman. Edge & Johnstone I 1 2005). 
Furthermore, since the southern radio jet is, by its brightness, in- 
ferred to be pointing towards us, we can conclude that the disc is 
seen in projection with its northern portions closer to us, and that it 
is rotating in a counter-clockwise direction. 

On larger scales the kinematics of H2 1-0 S(l) 2.122 \im and 
[Fe II] 1 .644 urn is more complex. Where detected, both lines show 
a similar line-of-sight velocity field. An exception to this is the re- 
gion directly south and west of the blueshifted portion of the disc, 



where the H2 velocity field shows indications of redshifted veloci- 
ties. The H2 emission from this region is interpreted as a streamer 
of molecular gas, which is associated with the morphological fea- 
ture shown in Fig. |4p. The H2 velocity dispersion is characterized 
by a secondary peak at about 0.4 arcsec west of the nucleus in 
the upper right panel of Fig. [5] At this position the observed H2 
line profile splits into separate kinematic components and becomes 
double-peaked. Since the single-Gaussian fit is not an appropriate 
model for these complex line profiles, we show the velocity slices 
for the H2 emission in Fig. [7] The dominant streamer feature is 
seen as an extension to the south-west, which is associated with 
redshifted velocities. As argued in Section [4] we suggest that this 
represents an infalling stream of molecular gas. In this scenario, the 
gas stream enters the nuclear region on a retrograde orbit with re- 
spect to the sense of rotation of the molecular disc. This is likely to 
cause shocks and turbulence that lead to energy dissipation in the 
disc, which may stimulate more rapid gas accretion onto the nu- 
cleus of NGC 1275. In addition to this prominent streamer, Fig. [7] 
shows indications of other streamers which are likely to connect to 
the filamentary structure around NGC 1275 on larger scales. 

A common feature of the kinematic maps in both, H2 1-0 S(l) 
and [Fell] 1.644 [im emission, is the circular region of redshifted 
velocities at a distance of about 1.2 arcsec from the nucleus to the 
north-west. This region is characterized by small velocity disper- 
sions in both lines. With respect to the core, this region is found at 
a position angle of ~ 310° — 340°, which is close to the projected 
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Figure 2. X-band spectrum of the core of NGC 1275 integrated over a 1.0 arcsec-diameter circular aperture. The wavelength scale is shown as observed for 
NGC 1275 at z = 0.017559. Note the broad base of the He I and Br7 lines, which presumably originate in the inner disc (i.e. within the H2 emitting region). 



axis of the radio jet at PA=160° jPedlar et al. Il990t) . (i.e. PA=340° 
for the northern radio jet with respect to the core). Since the pro- 
jected location of this region roughly coincides with the jet axis, 
this feature may result from a jet-cloud interaction. 



6 IONIZATION STRUCTURE 

By comparing the H 2 1-0 S(l), [Fell] 1.644 11m Hel2.058 vim and 
Br7 emission, we find evidence of an ionization structure within the 
circum-nuclear disc. H2 1-0 S(l) emission traces the outer regions 
of the disc and is likely to form a toroid, rather than continuing 
all the way to the nucleus. We use the dust sublimation radius as 
an estimate for the innermost radius at which H2 molecules can 
survive without being dissociated by the nuclear radiation field. A 
simple estimate of that radius is given by the distance from the core 
at which the nuclear bolometric luminosity, Lboi, is sufficient for 
heating the dust to its sublimation temperature, T au b: 

Lboi = 47rila ub aT s t lb , (2) 

where a is the Stefan-Boltzmann constant. We use 4 x 10 44 erg s - 1 
as bolometric luminosity of NGC 1275 jLevinson et al. I [l995h . 
For a dust sublimat ion temperature of ~ 1500 K (see e.g. 
iNenkova et al. I l2008h . the sublimation radius is about 0.1 pc. A 
larger sublimation radius is possible, if the sublimation tempera- 
ture is lower: For example, 1 pc corresponds to a local tempera- 
ture of about 500 K. For the forbidden [Fell] 1.644 urn line, the 



critical density is of the order of 10 s cm -3 (see Fig. 112b. There- 
fore, the [Fe II] 1 .644 [im emission becomes less important relative 
to the permitted lines in the high-density regions closer to the nu- 
cleus. The He I 2.058 (xm and Br7 lines, however, are expected to 
be tracers of the (presumably photoionized) inner portions of the 
disc as well as, possibly, the outer rim of an inflated disc facing 
the AGN. Evidence for such a scenario can be found in Fig. [8] 
This figure shows a single-Gaussian fit to the He I 2.058 |i.m and 
Br7 emission lines in NGC 1275, in direct analogy to Figs [5] and 
[6] for the H2 1-0 S(l) and [Fell] 1.644 \im emission lines. The 
He I 2.058 |J.m and Br7 emission is slightly extended compared to 
the /("-band continuum emission of NGC 1275 (see also Fig.O, but 
more compact than the H2 1-0 S(l) and [Fe II] 1.644 |j.m emission. 
The orientation of the major kinematic axis of the velocity field in 
He I 2.058 |im and Br7 matches that of the disc in H2 1-0 S(l) and 
[Fe II] 1 .644 |j.m, suggesting that the He I 2.058 (xm and Br7 emis- 
sion are associated with the same physical structure as the latter two 
emission lines. In addition, like H2 1-0 S(l) and [Fell] 1.644 |j.m, 
He I 2.058 |J.m and Br7 show a peak in velocity dispersion that is 
centred on the core of NGC 1275. However, the velocity disper- 
sion of He I and Br7 in the core (~ 330 km s _1 ) is larger than 
that of H 2 1-0 S(l) (~ 240 km s" 1 ). The same result can be de- 
duced from the integrated nuclear spectrum in Fig. [2] where the 
He I 2.058 !im and Bi'7 lines are shown to be broader than the H2 
lines (see also FWHM values in Table [TJ. This larger nuclear ve- 
locity dispersion of the He I 2.058 |J.m and Br7 emission compared 
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Figure 5. Composite result of fitting single-Gaussian profiles to the H2 1-0 S(l) 2.122 |xm emission in each NGC 1275 spectrum. The Gaussian central 
velocity and the Gaussian dispersion in km s — 1 are shown in the upper left and right panels, respectively. The lower-left panel shows the continuum flux 
density in units of 10 — 15 erg s — 1 cm -2 A -1 . The lower-right panel shows the intensity of the Gaussian fit to the H2 line in units of 10 erg s" 1 cm -2 . 
The location of the continuum emission peak is marked by a black cross in each panel. North is up and east is to the left. The scale bar in the velocity map 
indicates 100 pc. 



to the H2 emission indicates that the He I and Br7 emission arises 
at least partially from the inner portions of the circum-nuclear disc. 
These inner parts of the disc are expected to be dominated by high 
gas velocities that manifest themselves in the form of broad line 
profiles as a consequence of PSF smearing. Fig. [6] shows that the 
[Fell] 1.644 \im emission in the core of NGC 1275 is also charac- 
terized by a larger velocity dispersion - i.e. a broader line profile 
in Fig.Q]and Table[T|- than the H2 emission. The velocity disper- 
sion displayed by the nuclear [Fe II] 1.644 |j.m emission is compa- 
rable to the value found for the He I 2.058 |J.m and Br7 emission. 
This suggests that the core of NGC 1275 on scales smaller than the 
NIFS resolution also contains [Fell] 1.644 urn-emitting gas with 



large velocities from ordered or turbulent motion, which may be 
associated with the disc rotation or with jet-gas interactions. 



7 EXCITATION MECHANISMS 

The excitation mechanism for the H2 and [Fe II] emission in the 
centre of NGC 1275 is analysed based on the line flux measure- 
ments given in Table [T] 



7.1 Molecular hydrogen 



Krabbe et al. 



d2000h 



Wilm arTEdge & Johnstonel 



and, 



mo re 
(2005) 



recently, 
and 
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Figure 6. Same as Fig.|5]for the [Fe II] 1.644 |xm emission in each NGC 1275 spectrum. 



iRodriguez-Ardila. Riffel & Pastorizal d2005h have analysed the H2 
excitation mechanis m in the nucleus of NGC 1275. All find it to 
be thermally excited, [wilman, Edge & Johnstone I J2005h find that 
the excitation temperature is ~ 1360 ± 50 K from low-excitation 
transitions, with a higher te mperature of 3100 ± 500 K applying to 
high-e xcitation lines, while iRodriguez-Ardila. Riffel & Pastorizal 
(2005) have used the ratio of the H2 lines to determine a vibration 
temperature of 2200 ± 100 K, and a rotational temperature of 
1300 ± 100 K. 

The relative H2 level populations for the central 1.0 arcsec - 
diameter aperture centred on the nucleus of NGC 1275 are shown 
in Fig. [9] together with the corresponding measurements for a 0.5- 
1.0 arcsec-diameter annulus. The level populations are calculated 
from H2 line fluxes measured from the integrated //-band and K- 
band spectra (Figs[T]and[2]l and listed in TableQ] At densities above 
~ 10 5 cm -3 , the H2 gas is expected to be collisionally excited and 
so will be in thermal equilibrium. The level populations are then 



defined by an excitation temperature, T cx , that will be similar to 
the local gas temperature and can rise to several thousand Kelvin 
through heating by shocks or UV and X-ray irradiation. The relative 
level populations are derived by the Boltzman equation such that 
log(FX/Ag) = — Sup/Tex + constant, where F is the line flux, A 
is the line wavelength, A is the spontaneous emission coefficient, g 
is the statistical weight of the upper level of the transition, and E up 
is the energy of the upper level of the transition. These are plotted 
in the population diagram in Fig. [9] versus the upper-level energy. 
Transitions for a gas in thermal equilibrium lie along a straight line 
in this diagram with a slope of — 1/T ex . 

Krabbe et al71 



As 



found 



by 



(2000) 



and 



IWilman. Edge & Johnstone! J2005h . H2 transitions from the 
core of NGC 1275 with E up < 10 000 K are characterized by 
an excitation temperature of T cx ~ 1400 K, while transitions 
from higher energy levels are characterized by a hotter excitation 
temperature. The dashed straight line in Fig. [9] corresponds to the 
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T cx = 1360 K value found bv lWilman. Edge & Johnstone I J2005h 
for low-energy transitions. Our data are well fit by this tempera- 
tu re, although our low energy data ar e not as extensive as those 
of IWilman. Edge & Johnstone] J2005h who measured several Q- 
branch transitions that lie redwards of our spectra. A least-squares 
fit to only higher energy transitions (E up > 10 000 K) gives an 
excitati on temperature of T ex ~ 4290 K (solid straight line in 
Fig. |9>. IWilman. Edge & Johnstone] d2005h derive a comparable 
value of 3100 K from their data. 

A significant contribution from fluorescent excitation of the 
H2 emission can be ruled out based on the flux ratio of the high- 
energy H2 2-1 S(l) line and the low-energy H2 1-0 S(l) line, shown 
in a pixel-by-pixel comparison in Fig.[l0] The ratio is spatially uni- 
form at a value o f about 0.08-0 .09. Following the diagnostic dia- 
gram in fig. la o f M oun this value overlaps with the ther- 
mal excitation models but is much lower than the typical value of 
0.5-0.6 predicted for fluorescent excitation. Furthermore, the lack 
of spatial variation in the line ratio over the nuclear region contra- 
dicts a scenario with fluorescent excitation. If the higher excitation 
temperature of the high-energy levels were due to irradiation of the 
gas by nuclear UV or X-ray photons, this effect would then be ex- 



pected to be stronger close to the nucleus. In general, no significant 
change is observed in the behaviour of the high-energy transitions 
in our data between the 1.0 arcsec-diameter circular aperture and a 
0.5-1.0 arcsec-diameter annulus. This can be seen in Fig. [9] where 
the low- and high-energy transitions corresponding to the annulus 
(red symbols) display the same slope as those corresponding to the 
circular aperture (black symbols). 

The likely scenario for the excitation of the H2 emission in 
the nuclear region of NGC 1275 is shocks. In Fig. QT| we com- 
pare the two straight lines corresponding to the excitation temper- 
atures of the low and high-energy transitions in our H2 data for 
NGC 1275 to H2 data for shocked g as in the Orion molecular out- 
flow, taken from lBrand et al.l dl988l) . The shocked H2 gas in Orion 
shows the same properties as our H2 data for NGC 1275 in terms of 
a higher excitation temperature for higher-energy transitions, pre- 
sumably arising in post-shock gas having a range of temperatures. 
A range of excitation temperatures with higher temperatures for 
higher-energy transitions is also characteristic of the mid-infrared 
H2 excitation di agrams of NGC 1 275 and other H2-luminous radio 
galaxies, which lOgle et al. I feOlOh find to be consistent with shock- 
excitation. Therefore, we conclude that shocks are the dominant 
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Figure 8. Same as Fig.[5]for the He I 2.058 nm emission (top row) and the Bry 2. 166 )^m emission (bottom row) in each NGC 1275 spectrum. 
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Figure 4. (a) Integrated surface brightness of the H2 1-0 S(l) emission, 
(b) unsharp mask of the integrated H2 1-0 S(l) surface brightness, (c) inte- 
grated surface brightness of the [Fe II] 1.644 )^m emission, and (d) unsharp 
mask of the integrated [Fell] 1.644 |xm surface brightness. The images are 
plotted using an inverted grey-scale. The location of the respective con- 
tinuum peak in the K- and H-b&nd is indicated with a white cross. The 
H2 1-0 S(l) and [Fell] 1.644 u.m emission lines are integrated over ve- 
locity intervals from —500 km s" 1 to +500 km s — 1 and —800 km s — 1 
to +800 km s , respectively. The halo of negative pixels around the line 
emission is an artefact of the unsharp masking procedure, where a smoothed 
copy of the original image is subtracted from the original image (see text 
for details). The extended circum-nuclear emission in H2 is interpreted as 
a large-scale molecular accretion disc (panel a) on which a streamer is su- 
perimposed to the south-west of the nucleus (panel b). The streamer is not 
seen in the [Fe II] 1 .644 |im emission (panel d). 
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excitation mechanism for the H2 emission in the central 1.0 arcsec 
of NGC 1275. 



Figure 9. H2 population levels in the spectra integrated over a 1.0 arcsec- 
diameter aperture (black symbols) and a 0.5-1.0 arcsec-diameter annulus 
(red symbols) centred on the nucleus of NGC 1275. Relative upper-level 
population number densities are plotted versus the upper-level energy, ex- 
pressed as a temperature. H2 1-0 {circles), 2-1 {squares), and 3-2 {trian- 
gle) vibrational transitions are shown. Error bars are calculated for a con- 
stant flux error of ±2.5 X 10~ 16 erg s _1 cm -2 for all lines. The dashed 
straigh t line corresponds to the exci tation temperature of 1360 K derived by 
IWilman. Edge & Johnstone I fcOOSl) for transitions with E up < 10 000 K, 
normalized here to the H2 1-0 S(l) transition in the 1.0 arcsec-diameter 
aperture. The solid straight line shows a least-squares fit to the transitions 
with E up > 10 000 K from the 1.0 arcsec-diameter aperture. 



7.2 [Fe 11] 

Several authors have attempt ed to understand the nature of [Fe II] 
excitation in NGC 1275 (e.g. lRudv et al.|[T993h . The ionization po- 
tential of Fe + (16.2 eV) is similar to that of hydrogen, so [Fell] 
emission tends to arise in partially ionized regions. Being due to 
forbidden transitions between low-energy levels, the intensities of 
the near-infrared [Fe II] emission lines are principally dependent on 
the electron density. Calculations of the statistical equilibrium of 
collisional and radiative transitions among the lowest 16 levels of 
Fe II have been performed by several authors to study Seyfert galax- 
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Figure 12. Theoretical [Fe II] emission-line ratios versus electron density. The models are computed for a range of electron temperatures (1, 3, 5, 10, and 20 
X 10 3 K lowest to highest in plots). The shaded regions correspond to the measured line ratios and their uncertainties in the integrated //-band spectrum of 
NGC 1275 (Fig.[T). 



ies (e.g.. lMouri. K awara. & Taniguchi 2000; Thompson |l995h and 

!6J) 



Pcscnti et al 



2003 



Takami et al. 2006). The en- 



young stars (e.g., 
ergy level diagram for [Fe II] has been presented by these authors. 
Including only the lowest four terms is justified at the electron 
temperatures involved, which are likely to be < 20 000 K and 
leave higher energy levels significantly less populated. We have re- 
peat ed this 16-level atom calculatio n for Fe II using level energies 
from lAldenius & Johansson" I d2007l) . radiat ive transition pro babili- 
ties fr om lOuinet et al. I dl99~ ), (as li sted bv|Aldenius & Johansson I 



20071). and collision stren gths from lZhang & Pradhan 7T995h and 



Pradhan & Zhang I dl993l) . The model results for a range of [Fe II] 
emission-line combinations are shown in Fig.[T2]along with the em- 
pirical ratios and their uncertainties for NGC 1275 measured from 
our integrated nuclear spectrum (Fig. [TJ. The line ratios indicate a 
density in the [Fe II] region of ~ 4000 cm -3 . 

Deriving the electron temperature in the [Fe II] region is more 
difficult. The near-infrared [Fell] lines arise from transitions be- 
tween the lowest three terms (a 6 D, a 4 F, and a 4 D). The transitions 
detected in our //-band spectra all occur between the a 4 D and a 4 F 
terms and so have a limited range of upper-level energies. The J- 
band transitions of [Fell] (e.g., at 1.257 (j.m) share the same up- 
per energy levels, so do not provide temperature discrimination. 



Several authors have suggested using shorter wavelength transi- 
tions from the a 4 P term or longer wavelength transitions between 
the a 4 F and a 6 D terms to constrain the electron tem perature in 
Sevfe rt galaxies {Mouri. Kawara. & TaniguchTll2000|) and young 
stars dNisini et al.1 20021 ; Pesenti et alj|2003l) . This is the approach 
that we adopt here. 



We begin by using near-infrared transitions from the a 4 P 
energy level to constrain the electron temperature. The a 4 P 3 / 2 - 
a 4 D 7 / 2 transition occurs at 1.7489 p_m in the //-band so could 
make a nearly reddening-independent indicator with the 1.644 p.m 
line. The predicted dependence of this ratio as a function of electron 
temperature is shown in Fig.[T3] In practice, the [Fell] 1.7489 (im 
line coincides in wavelength with the H2 1-0 S(7) 1.7480 p.m line, 
which is strong in NGC 1275. Consequently, we have assumed that 
the [Fell] 1.7489 pm line is not detected in our //-band spectrum 
and have defined a 3a upper limit on the 1.644 p.m/1.749 p.m ra- 
tio accordingly. This is indicated in Fig. [T3] by the dashed line 
and solid arrow. This limit does not provide a useful constraint 
on the electron temperature. The [Fell] a 4 P3/ 2 -a 4 F 3 /2 0.9474 \im 
line, which has been sug gested as an elec tron temperature in- 
dicator for young stars dNisini et al. 1 12002|) . is similarly over- 
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Figure 10. Ratio of integrated H 2 2-1 S(l) 2.248 p flux to H 2 1-0 
S(l) 2.122 |xm flux. Each line is integrated over the velocity range ±250 
km s — 1 . An average continuum has been subtracted in each case. The loca- 
tion of the continuum emission peak is shown by a cross. North is up and 
east is to the left. The scale bar indicates 100 pc. 
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Figure 11. Comparison of the population diagram for H2 in the central 
1 arcsec of NGC 1275 with H2 data for shocked gas in the Orion molecu- 
lar outflow from Brand e t alj j 1 988t) - The straight lines are reproduced from 
Fig. [9] They correspond to the derived excitation temperatures in NGC 1275 
of 1360 K for transitions with E up < 10 000 K (dashed line) and 4290 K 
for transition s with E U y > 10 000 K (solid line). The data points are taken 
from lBrand e t al. 1 1988) and represent shocked H2 gas in the Orion molec- 
ular outflow. The symbols are the same as in Fig. [9] 
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Figure 13. Theoretical [Fell] emission-line ratios based on the a 4 P 3 / 2 - 
a 4 D 7 / 2 1-7489 p transition. Models are computed for a range of electron 
temperatures and electron densities, as indicated. The shaded region corre- 
sponds to the measured 1.644 u.m/1. 533 u.m ratio and its uncertainty in the 
integrated //-band spectrum of NGC 1275 (Fig. [TJ. The dashed line and 
arrow correspond to the limit placed on the 1 .644 u.m/1 .749 u.m ratio by the 
same data. 
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Figure 14. Theoretical [Fell] emission-line ratios based on the a 4 P 5 / 2 - 
a 4 Fg 12 0.8617 urn transition. Models are computed for a range of electron 
temperatures and electron densities, as indicated. The horizontal shaded 
region corresponds to the measured 1.644 u.m/1. 533 p ratio and its 
uncertainty in the integrated //-band spectrum of NGC 1275 (Fig. [TJ. 
The vertical dashed line and arro w correspond to the limit placed on 
the 1.2 57 u.m/0.862 u.m ratio from Riffel, Rodriguez- Ardila. & Pastoriza I 
120061) . 



whelmed by strong [S III] 0.9531 yim em ission in NGC 1275 
dRiffel. Rodn'guez-Ardila. & Pastoriza1l2006l) . 

iMo uri. Kawa ra. & Taniguchi ( 2000) for Seyfert galaxies and 
IPesenti et alj d2003l) for young stars have both suggested using 
the [Fell] a 4 P 5/ / 2 -a 4 F 9/ /2 0.8619 [im line strength as an indi- 
cator of electron temperature. We choose to compare this line 
to the [Fell] 1.257 [im line. This avoids reddening and aper- 
ture size issues to the greatest extent possible. Our model pre- 
dictions for the dependence of this line strength on electron tem- 
perature are shown in Fig. [14] The se results are in good agree- 
ment w ith the similar calculation of IMouri. Kawara. & Taniguchi I 
(2000) . We use the flux of the [Fell] 1 .257 \un line from 
Riffel, Rodriguez-Ardila. & Pastoriza 1 J2006h . A prominent emis- 
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sion line is also present at a rest-frame wavelength of ~ 
8630 A in the J/J-band spectrum o f NGC 1275 presented by 
Riffel, Rodr fguez-Ard iia. & Pastoriza I (2006). The feature is at a 
significantly different wavelength to the Ol 8446 A line seen in 
several others of their objects and to the closest hydrogen recom- 
bination lines (Hi P13 at 0.8667 and P14 at 0.8601 \im), and 
also the Call triplet emission at 0.8498, 0.8542, and 0.8662 ixm. 
We estimate the line flux to be ~ 2.0 ± 0.2 x 10" 14 erg s" 1 
cm -2 . If this is a detection of the [Fell] line, it indicates that the 
electron temperature in the [Fe Il]-emitting region is ~ 15 000 K 
(Fig. I14t . Otherwise, this represents an upper limit to the elec- 
tron temperature. [Fell] 0.8617 [im emission has been seen pre- 
viously in the Orion Nebula and in th e Seyfert galaxy NGC 4151 
l lOsterbrock. Shaw. & Veilleuxi ri990). which strengthens the case 
for a detection in NGC 1275. 

Possible excitation mechanisms for the [Fell] emis- 
sion in NGC 1275 are X-ray heating or fast shocks. 
iMouri. Kawara. & TaniguchTl ([2000) argue that [Fell] excita- 
tion by X-ray heating (i.e., a power-law photoionizing continuum) 
leads to higher electron temperatures (~ 8000 K) than shock 
heating (~ 6000 K). To the extent that their models are appli- 
cable, the elevated electron temperature in NGC 1275 is more 
indicative of X-ray heating than shock excitation. The high 
[Fell] 1.644 nm/Hl Bpy 2.166 ^m ratio of ~ 7.9 in NGC 1275 
also supports X-ray heating as the dominant excitation mechanism 
in the nuclear spectrum integrated over a 0.5 arcsec aperture. This 
ratio is a measure of the relative volumes of the partially-ionized 
region emitting [Fe II] and the fully-ionized region emittin g H I. It 



is exp ected to be in the range 0.1-1.4 in starburst regions l Colinal 
1993), but can have values up to ~ 20 i n environments where 
X-ray s penetrate deeply into neutral gas dAlonso-Herrero et al. I 



ray : 

1 19971) . The high value of this ratio in NGC 1275 and the proximity 
of the [Fe II] emission to the nucleus suggest that this is the process 
that excites the [Fell] emission in NGC 1275. 

We can also use th e J-band data of 

iRiffel. Rodrfguez-Ardiia, & Pastoriza] j2006h to locate NGC 
1275 in the e([FeIl] 1.257 nm)/e([Sll] 1.03 ixm) versus 
e([Fei i] 1.533 um)/e([Fen] 1.644 ixm) diagram of iPesenti et al.1 
(2003). They propose this as a diagnostic of Fe depletion 
on dust grains in the jets of young stars. Accord ing to ta- 
ble 5 in Riffel, Rodrfguez-Ardiia, & Pastoriza (2006]), the ratio 
[Fell] 1.257 nm/[Sll] (1.029 + 1.032 + 1.034 + 1.037) Jim 
in NGC 1275 is 0.5. In contrast, Galactic Herbig-Haro objects 
show a ratio of 3.0 ( IPesenti et al.l 120031) . These data, therefore, 
suggest that Fe is far more depleted in NGC 1275 than in Galactic 
Herbig-Haro objects. This is consistent with other estim ates for 
Seyfert galaxies (e.g., M ouri, Kawara, & Taniguchi 1 l200Ch where 
the strength of the [Fell] emission is interpreted to be primarily 
due to the ionization conditions, rather than to significant dust 
grain destruction in shocks. 



8 DISCUSSION 

8.1 Accretion scenario for the molecular disc 

The NIFS data for the circum-nuclear region of NGC 1275 
suggest predominantly shock-excited molecular-hydrogen emis- 
sion. This emission traces a rot ating disc, which, according to 
Wilman, Edge & Johnstone ( 2005), is likely to be unstable or to 
have a clumpy structure, as well as one or several streamer-like 
features. A possible interpretation of these results is an accretion 



scenario in which gas loses orbital energy as a result of shocks and 
turbulent dissipation, leading to a net mass inflow and to a net in- 
crease in binding energy of the accretion flow. Our model for this 
scenario is based on the model suggested for the ionized accretion 
disc in the centre of M87 by Do pita et al. I i 19971 see their equations 
(7) and (8)). In this model, at each radius in the accretion flow, the 
net increase in binding energy is matched by the shock luminosity 
from radiative losses, L(R): 



L(R) 



_d_ 

Jr. 



M(R)<t>(R) 



(3) 



where cj>(R) is the gravitational potential and M(R) is the net mass 
accretion rate including the inflow rate minus any mass outflow in 
a disc wind. If the fraction of the total shock luminosity that is 
radiated into the H2 line is given by k = St /Sh 2 > then 



L(R) = 2ttRkSh 2 



(4) 



First, we consider the case that the gravitational potential in 
equation (0 is dominated entirely by the central massive object 
(4>(R) = <j>Bti(R) = GMbh/R)- In this case, the net gain in 
binding energy in the accretion flow strongly depends on the radius 
and increases steeply towards small radii. Combining equations {3} 
and (0, results in 



2-7T R k Sh 2 , 



(5) 



MGMbh 
2R? 

if we assume that the net mass accretion rate M(R) = M is ap- 
proximately constant (as is likely in the outer parts of the accretion 
disc) and also that the ratio k for radiation into the line is approx- 
imately constant. This means that the radial variation in surface 
brightness should follow a steep scaling law with radius, propor- 
tional to R~ 3 : 



Sh 2 oc MMbhR" 



(6) 



The assumption of constant accretion rate and radiative efficiency 
per H2 molecule is clearly not valid for the innermost parts of the 
accretion disc. In these regions, presumably, disc winds lower the 
net rate of mass accretion, H2 molecules are dissociated, and the 
critical density may be exceeded. All these effects lead to a flatten- 
ing of the H2 surface-brightness profile at small radii. This can be 
expressed in the form of a core radius so that 

So 



Sn 2 (R) = 



(7) 



l + (iJ/a H2 ) 3 ' 

where So is the amplitude and an 2 is a core radius. 

Now, we consider the case that the gravitational poten- 
tial in equation ([3} is dominated entirely by an extended stel- 
lar mass distribution. In this case, the enclosed mass at each ra- 
dius increases with radius. The corresponding gravitational po- 
tential will lead to an accretion-powered H2 surface-brightness 
profile that is flatter than the one in equation JSJ. As the AGN 
dominates the near-infrared continuum emission of NGC 1275 
in the circum-nuclear region, the detailed stellar distribution 
within the PSF of the AGN is uncertain. However, in order 
to illustrate the accretion-flow scenario in a gravitational po- 
tential that is purely dominated by the stellar mass distribu- 
tion, we refer to the stell a r mas s- distribution model used by 
IWilma n. Edge & Johnstone I d2005l) . [wilman, Edge & Johnstone I 
d2005h fitted the stellar surface-brightness profile, obtained from an 
HSTNICMOS F160W image of NGC 1275, at radii > 0.5 arcsec, 
excluding the innermost AGN-dominated radii. They find that the 
stellar surface-brightness profile can be modelled with a modified 
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Figure 15. Fit to the surface-brightness profile of the integrated H2 1-0 S(l) 
emission of NGC 1275 shown in Fig. [4^. The surface-brightness profile 
(small plus signs) is obtained over an azimuthal range of 180°, avoiding 
the perturbed region south-west of the nucleus (see Fig.|4ji and b). The large 
diamonds show the linear interpolation of the profile onto a regular grid in 
radius. The solid red line shows the fit to the H2 profile using equation 1101 
convolved with the PSF model shown in Fig. [5] The fit results in a = 2.5 
and ajj 2 =0.17 arcsec. 



Hubble profile with a core radius of a, — 2.3 arcsec. For a con- 
stant mass-to-light ratio, the modified Hubble profile corresponds 
to a deprojected stellar density of 

-3/2 



p*(R) = pa 



1 



(£) 



(8) 



(cf. IWilman. Edge & Johnstone! l2005h . Our NIFS observations 
probe radii which are small compared to the above core radius of 
a, = 2.3 arcsec. On these scales, the stellar density from equa- 
tion l[8]l is approximately constant with radius (p*(J?) ~ po). This 
means that the enclosed stellar mass increases proportional to R 3 
and that the gravitational potential &*(R) is proportional to R 2 . 
According to equations ((3) and l(4j, this type of gravitational poten- 
tial results in an accretion-powered H2 surface-brightness profile 
that is constant with radius: 



Sh 2 (R) ~ constant. 



(9) 



In order to compare the above accretion-flow model for the 
H2 emission with the data for NGC 1275, we fit the observed H2 
surface-brightness profile with a generalized form of equation {7J, 
in which the power-law index a is included as a fit parameter: 

So 



Sh 2 (R) = 



(10) 



l + {R/a H2 ) a ' 

For a — 3, this function is equal to the model for the H2 surface- 
brightness profile in a gravitational potential that is entirely dom- 
inated by the central massive object (equation Q). For a = 0, 
this function is equal to equation (9}, which shows the example of 
an accretion flow in the gravitational potential of an extended stel- 
lar mass distribution. In order to take into account PSF smearing, 
equation J 10b is convolved with the PSF model shown in Fig. [3] 
before it is fitted to the observed data. The observed H2 surface- 
brightness profile for the inner 1 arcsec (358 pc) of NGC 1275 is 
displayed in Fig. Q3] The profile (small plus signs) is measured 
over an azimuthal range of 180° (in order to avoid the streamer 
to the south-west of the nucleus) and mapped onto a regular grid 
in radius via linear interpolation (large diamonds). If the incli- 



nation of the H2 disc caused significant projection effects, i.e. 
if the ellipticity of the isophotes were significant, the measure- 
ments for each surface brightness level would be distributed over 
a range of projected radii. As such projection effects are negligi- 
ble at log(i?) < —0.2 arcsec (225 pc) in Fig. 1 151 (small crosses), 
we do not deproject the profile. More significant projection effects 
are only seen at \og(R) > —0.2 arcsec. But the H2 distribution on 
these scales is likely to be more complex than the disc configuration 
inferred from the rotation pattern for R < 50 pc. 

The fit of the PSF-convolved function of equation J 10b to the 
observed H2 surface-brightness profile results in an 2 = 0.17 arc- 
sec and q = 2.5 (red line in Fig. 1 1 5b- This profile is slightly flatter 
than the profile expected for an accretion flow in a gravitational po- 
tential dominated by a central point mass. We, therefore, conclude 
that the observed H2 surface-brightness profile in NGC 1275 is in 
general agreement with the above accretion-flow model when tak- 
ing into account that the gravitational potential has contributions 
from an extended mass distribution in the inner 1 arcsec (358 pc). 

Wilman, Edge & Johnstone ] j2005h estimated that the en- 
closed stellar mass in NGC 1275 is about 1 per cent of the total 
dynamical mass at R = 50 pc (0.14 arcsec). It is therefore negli- 
gible on the scales that are characterized by the regular disc rota- 
tion (see Section [5}. As the surface-bright ness profile is fitted up 
to larg er scales, we use the model of I Wilman. Edge & Johnstone I 
J2005h in order to deduce that the enclosed stellar mass increases 
by a factor of ten from R = 50 pc (0.14 arcsec) to R = 110 pc 
(0.3 arcsec) and by a factor of 100 from R — 50 pc (0.14 arcsec) 
to R = 240 pc (0.67 arcsec). 

In addition to the stellar mass, the molecular gas mass is 
likely to contribute significantly to the total mass in the centre 
of NGC 1275. We derive a tentative estimate for the amount of 
molecular gas in the central R ~ 50 pc of NGC 1275 by com- 
paring our N IFS data to the l arge -scale H2 1-0 S(l) and CO ob- 
servations bv lLim et all d2012l) and lSalome e"t~ai] d2006|) .l Lim ct al. 
d2012l) report a total H 2 1-0 S(l) flux of the NGC 1275 nebula (nu- 
cleus plus filaments) of 3.87 x 10 -13 erg s _1 cm -2 . According 
to our NIFS data, of the order of 1 per cent of this flux originate 
from the central R < 50 pc. Assuming that this flux ratio translates 
into a mass ratio and that the same ratio applies to the near-, mid- 
infrared, and cold molecular gas, we estimate the cold molecular 
gas mass in the central R < 50 pc to be about 1 per cent of the total 
cold molecular gas mass of the NGC 1275 nebula (4 x 10 10 M , 
ISalome et ai]|2006t) . This results in a cold molecular gas mass of 
Mh 2 (R < 50 pc) ~ 4 x 10 s Mq, which is about two orders of 
magnitude larger than the stellar mass estimate for the same region. 

8.2 Disc model and black-hole mass 

As the H2 emission in the inner R ~ 50 pc of NGC 1275 shows 
evidence for a rotating disc, the enclosed mass of the disc can be 
inferred from kinematic modelling. If the potential of the rotating 
disc is dominated by the central supermassive black hole, the en- 
closed mass provides a direct estimate of the black-hole mass. For 



was reported by Wilman, Edge & Johnstone 1 fc005h . using the ob- 
served jump in H2 velocity across the nucleus as an indication of 
disc rotation. Their data, based on seeing-limited UIST (UKIRT 
Imaging Spectrometer) observations, lack the resolution to recon- 
struct the full two-dimensional intensity and velocity profiles of the 
disc. Their black-hole mass estimate is, therefore, based on a single 
value for the nuclear velocity discontinuity of 240 km s _1 resolved 
on a radial scale of 50 pc (~0. 15 arcsec). Another estimate for the 
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black-hole mass of NGC 1275 was given bv lBettoni et al.l J2003b . 
Using the M-a relation and a velocity dispersion of 250 km s _1 , 
they obtain an indirect mass estimate of 4.1 x 10 8 Mq. 

We have modelled the H2 emission for the inner 7? ~ 50 pc of 
NGC 1275 three-dimensionally, using a simulated NIFS data cube 
of a thin, transparent disc of H2 in Keplerian rotation around a cen- 
tral point mass. The model disc is first created in polar coordinates 
by assigning each radial and azimuthal sampling point a Gaussian 
line component. The central wavelength of the line component is 
given by the circular velocity at that radius, as determined from the 
Keplerian rotation curve. The line width is set to the NIFS spectral 
resolution of 60 km s _1 , i.e. the line components are assumed to be 
intrinsically unresolved. As discussed in Section [8T1 the observed 
H2 surface-brightness profile does not show any significant projec- 
tion effects on R ~ 50 pc scales. Therefore, we define the line 
flux at each radius in the face-on disc using equation J 10b with the 
parameters determined from the fit in Fig. [15] The resulting model 
disc is projected spatially onto an oversampled Cartesian grid of 
pixels. This step involves summing up all spectra that fall into each 
pixel for the given input inclination and the input position angle of 
PA = 68° (see Section [5}. Finally, in order to simulate the NIFS 
data, each two-dimensional slice of the model data cube is con- 
volved spatially, using the PSF from Section[3] and rebinned onto 
the rectangular NIFS pixel scale. The output is a full 3x3 arcsec 2 
NIFS data cube of an "observed" model disc. 

For the actual comparison between the model and the data 
via x 2 values, we only consider spectral pixels from the central 
~ 0.3 x 0.3 arcsec 2 (~ 110 x 110 pc 2 ) around the position of 
the nucleus of NGC 1275 (on which the model disc is centred). 
This ~ 110 x 110 pc 2 section covers the region with the most 
regular disc-type velocity field in the observations (see Fig. [5}. The 
surface-brightness profile (Fig.ll5t does not show a sharp cut-off 
at the edges of the kinematically defined disc but rather a smooth 
decline in H2 surface brightness over the full NIFS field-of-view. 
Therefore, we use an intrinsic size of the model disc that is much 
larger than the ~ 110 x 110 pc 2 region and extends well beyond 
this region for all projection angles considered here. For numerical 
purposes, the model disc is also truncated at small radii in order to 
avoid infinite velocities when the radius in the Keplerian disc ap- 
proaches zero. We choose an inner cut-off radius of 1 pc. This cor- 
responds to the minimum radius at which H2 molecules are likely 
to survive the nuclear radiation field, as outlined in Section[6] 

Fig.[T6]shows the \ 2 values for a comparison between the ob- 
served H2 1-0 S(l) lines and a grid of models with varying central 
mass and intrinsic disc inclination in the central ~ 0.3 x 0.3 arcsec 2 
(~ 110 x 110 pc 2 ) around the core of NGC 1275. In this figure, 
the absolute flux scale of each model has been adjusted indepen- 
dently to the value that fits the observed H2 lines in the central 
~ 0.3 x 0.3 arcsec best. As an example, the line profiles from 
one of the best models (i = 45° and M = 8.32 x 10 s M Q ) are 
compared to the observed line profiles in Fig.QT] It is obvious that 
the observed line profiles in the spectra to the west of the nucleus 
show a redshifted wing, which is not present in the models for pure 
disc rotation. This redshifted wing is part of the streamer discussed 
in Sections|4]and|5] 

Fig.[j6]shows that the enclosed mass of the best-fitting model 
depends on the disc inclination. This is due to the degeneracy be- 
tween central mass and inclination in reproducing the same line- 
of-sight velocities in the Keplerian disc. Small \ 2 values populate 
a valley extending from a central mass of about 2.1 x 10 9 M 
at an inclination of 25° to a central mass of about 7 x 10 8 Mq 
at an inclination of 65°. The global lower mass limit for the cen- 



3. 021 
2.5ll 
2.09I 
I.74I 
I.45I 
1.20 
1.001 
0.83 
0.69 
0.58 
0.48 
0.40 



■ 200 



I 



180 



25.0 



160 



140 . 



120 ■ 



100 



80 



60 



35.0 45.0 55.0 
inclination [deg] 



65.0 



Figure 16. Comparison between the disc models and the observed H2 data 
for a grid of models with varying central mass and inclination. The re- 
duced x 2 values are calculated based on a three-dimensional comparison 
between the modelled and observed data cubes. Spectrally, the x 2 is com- 
puted over a velocity range from —500 km s _1 to +500 km s _1 around 
the rest wavelength of H2 in NGC 1275. Spatially, the x 2 values are com- 
puted for boxes of 3 x 7 non-square NIFS pixels (i.e. ~ 0.3 X 0.3 arcsec 2 
or ~ 110 X 110 pc 2 ) centred on the nucleus, as displayed in Fig. 1171 The 
white contour level shows A\ 2 = 9 with respect to the minimum x 2 - The 
black vertical lines mark the inclination range of 45° ± 10° assumed for 
the black-hole mass estimate. 



tral mass of > 5 X 10 8 M in Fig.[l6]is robust against a possible 
underestimate of the PSF width in our K-band observations. By 
increasing the FWHM of the model PSF, the best-fitting models at 
a fixed inclination are shifted toward higher central masses. This is 
a consequence of enhanced PSF smearing which causes the central 
velocity gradient in the models to decrease, so that a higher central 
mas s is required in order to fit the observ ed velocities. 

Wilman, E dge & Johnstonel d2005h use an inclination of 45° 
in their estimate for the dynamical mass in NGC 1275. This value 
is based on the assumption that the disc rotation axis has the same 
orientation as the axis of the radio je t, which is inclined by 30°- 
55° w ith respect to the line-of-sight (Walker, R omnev & Benson] 
1 1994b . For a similar disc inclination of 45° ± 10°, our model sug- 
gests a best-fitting central mass of M con = (8 ±1) x 10 8 M 
within the contour level sho wn in F ig. [T6l T his is larger tha n 
the 3.4 x 10 8 M Q derived bv lWilman. Edge & Johnstone I J2005h . 
This discrepancy is a consequence of mainly two factors: First, 
the higher resolution NIFS data resolve a larger velocity discon- 
tinuity across the nucleus. Secondly, in contrast to the calcula- 
tion by 1 Wilman. Edge & Johnstone l( l2005h . our three-dimensional 
analysis takes into account PSF smearing, which leads to a 
decrease in the observed velocities close to the nucleus. Our 
value for the central mass is in agreement with the M-a rela- 
tion within the rms scatter. Assuming a stellar v elocity disper- 
sion of a = 250 kms -1 for NGC 1275 (see iBettoni et al. I 
120031). the updated M -a relation for all morphological types from 
lGrahametal.H201lb .log (M BH /M ) = (8.13±0.05) + (5.13± 
0.34) log [a/200 km s" 1 ] results in Mbh = 4.2 x 10 s Mq, which 
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Figure 17. Comparison between the H2 profiles in the data (black) and in the disc model for i = 45° and M = 8.32 X 10 Mq (red). The figure shows the 
central 3x7 pixels (i.e. ~ 0.3 X 0.3 arcsec 2 or ~ 110 X 110 pc 2 ). In each spectrum, the dotted vertical line indicates the systemic velocity of NGC 1275. 
Tick marks on the spectral axis are shown in steps of 200 km s^ 1 . The non-square spatial pixels of NIFS have scales of 0.044 arcsec per pixel in x and 
0.103 arcsec per pixel in y. Observational artefacts affecting some of the observed spectra are excluded from the \ 2 calculation in Fig.[l6]and marked by a 
light gray colour. North is up, east is to the left. 



is consist ent with our value wi thin the rms scatter of 0.43 dex re- 
port ed bvlGraham et al.l d201 ll). 

IWilman. Edge & Johnsto ne (2005) show that the stellar con- 
tribution to the total dynamical mass within R ~ 50 pc is negli- 
gible and, therefore, argue that the dynamical mass represents the 
black-hole mass. However, as shown in Section [8T1 we estimate 
that a large reservoir of molecular gas of the order of Mh 2 ( R < 
50 pc) ~ 4 x 10 s M Q is present in the central parts of NGC 1275. 
This large gas mass could potentially contribute significantly to the 
derived enclosed mass. In fact, if this gas mass is subtracted from 
the total enclosed mass of (8 ±1) x 10 8 M Q , the resulting black- 
hole mass of ~ 4 x 10 8 Mq is almost exactly the same as the one 
derived from the M-a relation. However, as this result involves 
large uncertainties, we conclude by merely noting that our estimate 
of the enclosed mass may rather represent an upper limit for the 
black-hole mass. 



8.3 Comparison to other H2 -luminous radio galaxies 

Using Spitzer observations of p ure rotational H2 transitions in the 
mid-infrared, lOgle et al] detect warm molecular hydrogen 



emission in 30 per cent of their sample of z < 0.2 3C radio 
galaxies, including NGC 1275. These molecular hydrogen emis- 
sion galaxies (MOHEGs) are found among radio galaxies of both 
Fanaroff-Rile y classes and amo ng cool core clusters. The sam- 
ple studied bv lOgle et al.l ([20 10) shows that radio MOHEGs as a 
class typically have low star-formation rates despite the large H2 
masses, which indicates that turbulence prevents the gas from ef- 
ficiently forming stars. Furthermore, most MOHEGs harbour low- 
luminosity AGN that cannot provide sufficient X-ray heating for 
powering the observed H2 emission. The AGN in NGC 1275 
(logLx(2 - 10 keV) = 42.91) belongs to the more X-ray lu- 
m inous AGN amon g the radio MOHEGs (cf. fig. 17 and table 9 
in lQgleetalj[2O10l) . This means that X-ray heating may make a 
stronger contribution in NG C 1275 than in many X-ray faint MO- 
HEGs. However, fig. 17 in lOgle et alj jioich shows that the ra- 
tio of the observed warm H2 luminosity to X-ray luminosity for 
NGC 1275, i.e. L(R 2 )/L X = 0.07, is still a factor of 7 higher 
than the maximum ratio of 0.01 derived for an X-ray dissociation 
region. Therefore, even in NGC 1275, X-ray heating by the AGN 
is unlikely to be the primary excitation mechanism for H2. 

Instead, lOgle et alj JioToh conclude that the most likely exci- 
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tation mechanism for the H2 emission in radio MOHEGs is shocks 
(or, less likely, cosmic ray heating). The same mechanism is sug- 
gested for NGC 1275, based on the NIFS data discussed in this 
paper (Section [7.U . This raises the question whether these shocks 
originate from radio-jet fee dback, as su ggested for the overall sam- 
ple of radio MOHEGs bv lOgleetal] (2010), or from gas accre- 
tion, as discussed for the molecular disc in NGC 1275 in Sec- 
tion 18.11 In the feedback scenario, shocks are driven into dense 
clouds in a porous interstellar medium by the hot cocoon pow- 
ered by the radio je t , as fo r exa mple shown in the sim u lation s by 
IWagner & Bicknelll d201ll) and ISutherland & Bicknelj d2007l) . In 
the accretion scenario, gas clumps are repeatedly shocked while 
falling into the centre of the galaxy so that gravitational potential 
energy is converted into shocks and turbulent heating. In a case 
study of 3C 326 N, a z = 0.09 radio M OHEG with a very low 
star-formation rate and an X-ray faint AGN. IOgle et al.l d2007l) sug- 
gest a tidally-i nduced accretion flow a s the most likely origin of the 
shocks, while iNesvadba et al.l d2010h favour the jet feedback sce- 
nario. 

Simple energy considerations for NGC 1275 show that ei- 
ther of the processes, radio jet feedback or gas accretion, is suf- 
ficient to power t he observed mid- or near-infrared H2 emission. 
For the radio jet, lOgle et al] d20ld) report a jet cavity power of 
Pint Cavity = 1.4 x 10 44 ergs -1 , which is almost 3 orders of 
magnitude larger tha n the luminosity of the warm H2 emission 
(10 41 ' 75 erg s _1 , cf. lOgle et al.ll201Ct) . For the accretion process, 
we estimate an accretion luminosity of the order of 10 43 erg s" 1 , 
which is about one order of magnitude larger than the warm H2 
luminosity. This accretion luminosity is derived by dividing the po- 
tential energy of the infalling gas of about 3 x 10 58 erg by the 
dynam ical free-fall time of ~ 1 x 10 s yrs, in analogy to lOgle et al.1 
d2007l) . For this calculation we assume that the total cold gas mass 
of 4 x 10 10 M dSalome et al]|2006l) falls back from an average 
distance of 30 kpc into the potential well corresponding to a stel- 



the MBH-M bulge / host relations at 2 = by dBennert et alJl20lH) . 
using the M ccn from Section l8^2l as the black-hole mass. 

As neither jet feedback nor gas accretion stands out as the 
dominant energy source in NGC 1275, it is likely that the observed 
properties of gas excitation and kinematics are a consequence of 
both processes. Rotating molecular discs, as found in the circum- 
nuclear region of NG C 1275, or on scales of 3 kpc in 3C 326 N 
dNesvadbaetalj201ll) . suggest that accretion processes play a role. 
Since the radio jet cannot impart angular momentum on a gaseous 
disc, disc rotation is a sign of gas infall. This gas inflow may be 
associated with tidal galaxy-galaxy interactions, or with accretion 
from the h alo. These processes have been discuss ed for the case of 
3C 326 N dQgle et alj2007l ; lNesvadba et alj2010h . For NGC 1275, 
it is most likely that the gas inflows consist of ga s which falls back 
from the radio jet bubbles dSalome et af] |2008a). Such back flows 
of gas uplifted during the active jet phase are also seen in simu- 
lations (Sutherland & Bicknell 2007). Tidally-induced gas inflows 
could result from an interaction with the high- velocity system. But 
this scenario is less likely, since t he high- velocity syst em is still at 
a large distance from NGC 1275 dGillmon et a"il2004l) . 

Observational evidence suggests that turbulence is a com- 
mon c haracteristic of the mol ecular gas component in radio MO- 
HEGs. |Nesv^aeta[| d201lh found the molecular gas in the disc 
of 3C 326 N to be highly turbulent with line widths of FWHM ~ 
500 kms -1 . Broad H2 lines are also characteristic of the sam- 
ple of_jadio_galaxieswith fast ionized/atomic outflows, studied 
by iGuillard etail d2012h . Integrated over the full field-of-view, 



our NIFS data of NGC 1275 show an H 2 1-0 S(l) line width of 
FWHM ~ 350 km s -1 , which is comparable to the above val- 
ues. A fraction of this line width is likely to be attributed to or- 
dered streaming motions, which are characterized by small (< 
100 km s _1 ) and largely positive velocities (see Fig.UJ. Since the 
models for the circum-nuclear disc in NGC 1275 (Section [8T2l pro- 
vide a good representation of the data without including any signif- 
icant intrinsic line width, the velocity dispersion of the H2 1-0 S(l) 
lines are likely to be < 100 km s _1 (i.e. FWHM < 240 km s" 1 ). 
(The total widths of the H2 lines in the small-aperture spectra, listed 
in TableQ] are much larger than that, but are dominated by the disc 
rotation.) The fact that the dispersion-related line width in the disc 
is smaller than the total integrated line width in the NIFS field- 
of-view suggests that the accretion flow is energetically capable of 
producing the turbulent motions in the disc. Alternatively, the tur- 
bulence may be driven by the radio jet and its hot cocoon. Turbulent 
heating via radio jet feedback can affect a large volume if the gas 
is distributed in the form o f clumps in a porous interstellar medium 
dWagner & Bicknelll201 lh . 

Turbulent heating of the molecular gas component by the jet 
interaction is only compatible with the observed disc configurations 
if the jet is, at the same time, inefficient in expelling the molecular 
gas from t he disc. Such a pictur e is also suggested by the obser- 
vations by Guillardet al. (2012), who find that in radio galaxies 
with fast ionized and atomic outflows the bulk of the mid-infrared 
H2 emission does not participate in the ou tflow, despite showing 
signs of turbulent heating. Simulations by Suth erland & Bicknelll 
(2007), which probe the effects of a jet interaction with an almost 
Keplerian and turbulently supported disc of a warm clumpy inter- 
stellar medium in an elliptical galaxy, show that the jet interaction 
is, indeed, inefficient in driving the gas out of the disc region. How- 
ever, there is a significant gas transport from the inner to the outer 
portions of the disc, which will eventually leave the inner portions 
depleted. If the disc is porous, the mass transfer rates are lower than 
in the case of a homogeneous disc of gas. 

We estimate the porosity of the circum-nuclear disc in 
NGC 1275 by comparing the observed mean density of H2 
molecules in the central R ~ 50 pc with the density expected for 
the corresponding gas component from the ideal equation of state. 
As shown in Section [8Tl we tentatively estimate the total molecular 
gas mass within a radius of R ~ 50 pc to be Mh 2 (R < 50 pc) ~ 
4x 10 s Mq. The corresponding mean number density of molecules 
in the central volume is nu 2 = Mu 2 (R < 50 pc)/(2m p V) ~ 
2 x 10 4 cm -3 , where m v is the proton mass and V is approx- 
imated by the spherical volume enclosed by a radius of 50 pc. 
This value is similar to the clump density derived from the ideal 
equation of state, if we assume that the clumps are composed of 
cold molecular gas at a mean temperature of T = 50 K and 
that they are in pressure equilibrium with the hot ISM pressure of 
p/k = 10 6 cm -3 K, following the scenario by Wagne r"& Bicknefl] 
d201ll) . The similarity between the mean observed density and the 
estimated clump density from the ideal equation of state suggests 
that the cold molecular gas, if distributed in clumps, is likely to 
have a high filling factor. An independent measure of the actual gas 
density is provided by detec tion experiments for high-density gas 
tracers. lSalomeet al. (2008b) report the detection of HCN(3-2) in 
the centre of NGC 1275, which indicates a gas density of at least 
10 4 cm -3 . This value agrees with the mean gas density estimated 
above. But, as a lower limit for the actual gas density, it does not 
rule out the possibility of smaller filling factors. Detections of other 
high-density gas tracers are needed in order to obtain a better esti- 
mate of the filling factor. 
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The relative importance of accretion and jet feedback in ex- 
citing and heating the molecular gas in NGC 1275 and other radio 
MOHEGs is poorly constrained with the current data. Improved 
observational parameters for the molecular gas phase - e.g., from 
ALMA observations of high-density tracers - together with refined 
simulations are required in order to explore the effects of the jet in- 
teraction with the interstellar medium in these molecular-hydrogen- 
rich radio galaxies. 



9 SUMMARY AND CONCLUSIONS 

The paper is based on near-infrared integral field spectroscopy of 
NGC 1275, obtained with the NIFS instrument at Gemini North at a 
spatial resolution of ~ 0.1 arcsec. We have presented a scenario in 
which the shock-excited molecular hydrogen, traced via the near- 
infrared ro- vibrational H2 lines, enters into the circum-nuclear re- 
gion in the form of one or more streamers and settles into a turbu- 
lent, clumpy or unstable accretion disc close to the nucleus, where 
energy dissipation via turbulence and shocks drives the further ac- 
cretion process towards the core. The observed disc kinematics has 
been used as a tool to infer a direct estimate for the central mass 
enclosed by the disc via model data cubes of discs in Keplerian ro- 
tation. We have compared the results to other H2 -luminous radio 
galaxies and discussed the relative importance of jet feedback and 
accretion as a mechanism for shock-excitation and heating of the 
molecular hydrogen. 

In summary, the main results presented in this paper are: 

• The circum-nuclear H2 1-0 S(l) and [Fe II] 1.644 um lines in 
NGC 1275 show evidence for a disc on R — 50 pc scales, both in 
morphology and kinematics. The projected major kinematic axis of 
this disc is perpendicular to the projected axis of the radio jets. 

• A possible ionization structure within the disc is discussed by 
comparing the H2 1-0 S(l) and [Fell] 1.644 um emission with the 
He I 2.058 um and Br7 lines. The latter two lines are only found 
in a compact region around the core, but their kinematics is char- 
acterized by the same disc signature as seen in H2 1-0 S(l) and 
[Fell] 1.644 am. It is suggested that both the He I 2.058 am and 
Br7 emission lines arise at least partially from the inner parts of the 
disc, because they display a larger central velocity dispersion than 
the H2 1-0 S(l) line. Future higher angular resolution observations 
are required to fully explore this ionization structure. 

• The H2 surface-brightness profile in the inner 1 arcsec around 
the core is well represented by a PSF-convolved profile following 
an R~ 2 5 power law with an inner core radius of 0.17 arcsec. This 
type of surface-brightness profile is found to be in general agree- 
ment with an accretion model in which the line emission matches 
the net increase in binding energy of gas falling into the combined 
gravitational potential of a central massive object and an extended 
mass distribution. 

• In H2, the data reveal an elongated perturbation to the south- 
west of the nucleus, which is characterized by redshifted velocities. 
This perturbation is interpreted as the strongest of several molecular 
gas streamers, which are likely to be falling into the core. 

• The relative line strengths of the ro-vibrational H2 lines in the 
central 1 arcsec of NGC 1275 indicate thermal excitation of the 
molecular hydrogen via shocks. In agreement with previous stud- 
ies, we find excitation temperatures of 1360 K and 4290 K for tran- 
sitions with E upp < 10 000 K and E upp > 10 000 K, respectively. 
Any significant contribution from fluorescent excitation can be ex- 
cluded based on the small ratio of the H2 2-1 S(l) 2.248 um to 
the H2 1-0 S(l) 2.122 um flux and its spatial uniformity over the 



nuclear region. Instead, the H2 excitation in NGC 1275 shows the 
same characteristics as shocked gas in the Orion molecular outflow. 

• The high electron temperature as well as the high 
[Fell] 1.644 um/Hl Br7 2.166 urn suggest that the [Fell] lines in 
the integrated nuclear spectrum are excited by X-ray heating. Al- 
ternatively, the line emission can be excited by fast shocks. Models 
for a range of [Fe II] emission-line ratios are presented. They indi- 
cate a density of ~ 4000 cm -3 in the [Fe II] emitting region. The 
electron temperature in this region is less well constraint. Using 
less certain detections of [Fe II] lines at shorter wavelengths from 
previous publications, an electron temperature of ~ 15 000 K is 
estimated. 

• An estimate for the mass enclosed by the circum-nuclear H2 
disc of NGC 1275 is derived by modelling the data via simulated 
NIFS data cubes. Assuming a disc inclination of 45° ± 10°, the 
enclosed mass is estimated to be (8 lj) x 10 8 Mq, Interpreted as 
a direct measurement of the black-hole mass, this value is larger 
than previous black-hole mass estimates, but is consistent with the 
M-cr relation within the rms scatter. A tentative estimate of the 
molecular gas mass in the central region of NGC 1275 results in 
Mu 2 (R < 50 pc) ~ 4 x 10 s M , which is a non-negligible 
fraction of the enclosed mass. The enclosed mass may, therefore, 
represent an estimate for the upper limit of the black-hole mass. 

• Basic energy considerations show that the shock-excitation 
and turbulent heating of the molecular gas in NGC 1275 could be 
driven by the jet interaction, the accretion process, or a mixture 
of both. An order-of-magnitude estimate of the distribution of the 
cold molecular gas in the circum-nuclear disc of NGC 1275 indi- 
cates a high filling factor. Improved constraints on the properties of 
the molecular gas phase can be expected from future observations 
(e.g., ALMA). 
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